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Preface 

This  report  shews  the  results  of  my  experimental  determination  of 
the  actual  size  distribution  of  aerosols  collected  by  an  electrostatic 
precipitator  developed  at  the  Air  Force  Institute  of  Technology.  The 
specific  precipitator  investigated  is  a  model  which  has  evolved  over 
the  last  two  academic  years  and  which  has  been  used  for  the  collection 
of  particulate  matter  from  atmospheric  air  in  the  detection  of  its 
radioactivity  content.  The  ultimate  value  of  this  study  is  to  provide 
more  evidence  which  would  substantiate  continued  investment  of  time 
and  funds  to  perfect  an  electrostatic  precipitator  for  Air  Force  field 
use.  The  value  of  such  a  piece  of  equipment  goes  without  saying.  I 
feel  that  the  precipitator  investigated,  on  the  strength  of  Baker's 
work  (Ref  1)  and  the  data  presented  here,  warrants  continued  develop¬ 
ment  towards  a  working  model  for  field  applications. 

I  wish  to  thank  Captain  Charles  J.  Bridgman,  i^y  thesis  advisor, 
for  his  guidance  and  help  in  this  project.  His  enthusiasm  and 
interest  in  this  study  has  been  a  constant  source  of  encouragement 
to  me.  I  would  also  like  to  thank  Mr.  William  Schoonover  of  the 
Aerial  Reconnaissance  Laboratory  and  Major  William  Metscher  of  the 
Electronics  Technology  Laboratory  for  providing  the  necessary  space 
and  facilities.  I  am  grateful  to  Mr.  John  Blasingame  for  his 
generous  help  in  the  field  of  microscopy.  I  also  thank  CWO  James. 
Miskimen,  Staff  Sergeant  Floyd  and  Mr.  Elworth  of  the  Physics  Laboratory 
and  Mr.  Winston  Wolfe  of  the  AFIT  Workshop  for  their  assistance. 
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Abstract 

An  electrostatic  airborne  radioactivity  collector  (EARC-I)  is 
analyzed  bo  determine  what  correlation  exists  between  the  Banner  in  which 
it  actually  collects  atmospheric  aerosols  and  the  theoretical  model 
previously  established.  Size-distributions  of  aerosols  with  radii  in 
■the  range  of  1.0  to  10.0  microns  are  determined  by  optical  microscopic 
analysis  of  particulate  natter  collected  at  nine  locations  along  the 
collection  plates  of  the  precipitator.  Theoretical  size-distributions 
are  calculated  for  the  same  locations  based  on  the  IBM  7090  solution  of 
theoretical  equations  describing  the  charging  processes  in  EAIiC-I.  It 
is  concluded  that  i5AHC-I  does  collect  particles  roughly  in  accordance 
with  the  theoretical  model. 
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DETERMINATION  OF  THE  SIZE  DISTRIBUTION  OF  AEROSOLS  COLIECTED 
BY  ELECTROSTATIC  PRECIPITATION 

I,  Introduction 

The  purpose  of  this  study  is  to  determine  if  the  electrostatic 
precipitator,  EARC-I,  developed  at  the  Air  Force  Institute  of  Technology, 
does  in  fact  collect  atmospheric  aerosols  exactly  as  predicted  by  Lamberson 
(Ref  7)  or  nearly  so.  A  limited  portion  of  the  particle  size  spectrum 
is  considered.  Using  die  micron,  u,  as  the  unit  of  measurement,  the 
range  of  consideration  includes  those  particles  with  a  radius  equal  to 
1.0  u  up  through  particles  with  a  radius  equal  to  10.0  u.  The  terms; 
aerosols,  particles,  particulate  matter,  dust  particles  are  considered 
to  have  the  same  meaning  in  this  report. 

This  investigation  is  experimental  in  nature.  In  essence,  it  consists 
of  collecting  particulate  matter  from  the  ground- level  atmosphere  at 
Wright-Patterson,  AFB,  Ohio.  Samples  from  the  collected  matter  will  be 
analyzed  by  optical  microscope  to  determine  a  particle  size-distribution 
at  several  locations  along  the  collection  plates  of  EARC-I.  Using 
Lamberson' s  equations,  the  theoretical  size-distributions  will  be  determined 
for  the  same  locations.  An  attempt  will  be  made  to  explain  any  discrepancies 
between  the  theory  and  experimental  results.  Finally,  recommendations  for 
future  development  of  EARC-I  will  be  made  based  upon  the  experimental 
findings. 

The  value  of  this  investigation  is  in  such  recommendations.  If  the 
theory  and  experimental  data  agree  within  reason,  the  future  development 
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of  EAEC -I  is  justified  towards  the  eventual  use  by  the  Air  Farce  as  well 
as  by  Civil  Defense  agencies  far  radioactivity  sampling.  If,  on  the 
other  hand,  the  expert  mart  al  data  indicates  that  EARC-I  produces  adequate 
results  on  the  basis  of  unexplained  phenomena,  the  continued  investigation 
in  this  Held  may  not  be  justified. 

Background 

The  collection  and  analysis  of  radioactive  aerosols  is  a  project 
which  has  been  under  investigation  by  graduate  Nuclear  Engineering 
students  of  the  Air  Force  Institute  of  Technology.  A  composite  theoretical 
model  of  the  physical  processes  involved  in  electrostatic  collection  was 
developed  by  Lamberson  (Ref  7) .  A  high  volume  electrostatic  precipitator, 
EARC-I  (Electrostatic  Airborne  Radioactivity  Collector  Number  1),  was 
designed,  built  and  used  by  Baker  in  his  analysis  of  weapon  fallout  at 
Wright-Patterson  AFB,  Ohio  (Ref  l).  Stuart  (Ref  13)  did  a  design  optimization 
on  the  basis  of  the  theoretical  model. 

The  feasibility  of  electrostatic  precipitation  as  a  method  of  collecting 
radioactive  fallout  has  been  demonstrated  by  Baker.  In  gpite  of  this,  and 
the  long-standing  use  of  electrostatic  precipitators  for  various  industrial 
applications,  an  examination  of  the  literature  relating  to  this  subjec'. 
indicates  that  a  precise  knowledge  of  the  process  is  lacking.  Many  of  the 
processes  involved  are  at  best  described  only  after  a  host  of  simplifying 
assumptions  are  made.  As  a  step  towards  understanding  these  processes 
more  fully,  this  thesis  investigates  the  collection  of  naturally  occurring 
particles  in  such  a  precipitator. 
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II.  Apparatus 

EARO-I 

EARG-I  is  a  two-stage,  parallel-plate  electrostatic  precipitator 
which  is  the  actual  piece  of  equipment  desigied  and  constructed  by  Baker 
(Ref  1)  in  his  thesis  work.  The  only  work  done  on  the  precipitator  for 
■this  ejqperiment  was  the  replacement  of  all  resistors  in  the  current- 
limiting  network  and  the  rewiring  of  all  grids,  making  use  of  the  same 
materials  as  in  the  original.  So  there  were  no  design  changes  involved. 
The  exact  construction  details  sad  specifications  are  available  in  Bakers 
thesis  (Ref  1:3U).  Only  a  brief  description  and  sketch  are  included 
here  (Figure  2.1). 

In  the  form  of  a  rectangular  parallelopiped,  EARC-I  has  two  basic 
sections.  The  charging  section,  approximating  a  line-to-plate  geometry, 
is  a  wire-to-plane  arrangement .  The  plane  is  provided  by  l/8  inch  solid 
aluminum  sheeting.  The  wire  is  0,00$  inch  diameter  tungsten  wire  ti^itly 
drawn  across  an  aluminum  frame.  In  the  charging  section,  there  are  17 
layers  of  plate-grid -plate-etc. ,  alternating.  It  is  here  that  the  corona 
discharge  occurs  to  provide  the  mechanism  for  charging  aerosol  particles. 
Plate-grid  separation  is  2.0  cm.  Use  is  made  of  the  negative  corona, 
which  means  that  the  grid  (wire)  is  negative  with  respect  to  the  grounded 
plate. 

The  collection  section  consists  of  plate-to-plate  geometry  in  the 
form  of  17  horizontally  arranged  l/8  inch  solid  aluminum  sheeting  at 
2.0  cm  separation.  Nine  of  the  plates  in  the  collection  section  are 
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physical  extensions  of  the  plates  in  the  charging  section, 


The  microscope  used  is  a  Vicker ' s  Projection  Microscope  manufactured 
by  the  Research  Laboratories  of  Messrs.  Viekers-Arrastrongs,  Ltd.,  York, 
England.  (Figure  2.2) 

Lens.  The  following  optical  lens  were  used  throughout 

a.  Objectives  Cooke,  apochromatic,  oil-immersion  lenB,  with 
2.2  mm  focal  length,  numerical  aperture  of  1.32. 

b.  Projection:  Cooke,  1$  compensated. 

Lighting.  For  the  high  magnification  used,  advantage  was  taken  of 
the  carbon  arc  lanp  available  with  the  Vickers.  This  provided  a  high 
intensity,  evenly  distributed  circular  ligit  source.  The  lamp  con¬ 
denser  is  mounted  close  to  the  carbon  arc.  Adjustment  is  provided 
by  an  adjustable  iris  diaphragm  for  controlling  the  aperture  of  the 
lens.  The  lens  can  also  be  made  to  slide  toward  or  away  from  the 
light  source.  Figure  2,3  shows  the  image  forming  system.  The  working 
image  was  projected  on  the  ground  glass  projection  screen.  The  screen 
had  a  working  field  etched  with  a  10  cm  x  10  cm  area  which  was  further 
subdivided  into  1  cm  squares. 


5 


Figure  2.2 

Vicker's  Projection  Microscope 

1.  Ground-Glass  Projection  Screen 

2.  Carbon  Arc  Lamp 
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(Ref  17:  ft) 

Figure  2.3 


Image-Forming  S/stem  of  Vickers 
Projection  Microsoope 
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III.  Theory 

General 

The  theory  which  formed  the  basis  for  the  realization  of  EARC-I, 
although  presented  by  Lamberson  (Ref  ?) ,  is  in  general  agreement  with 
other  well-known  investigators  in  this  field,  e.g.  White,  Lowe,  Lucas, 
Penney.  The  chronological  sequence  of  processes  involved  are  briefly 
stated,  then  explained  in  more  detail.  Figure  3.1  should  be  an  aid  in 
visualizing  the  processes  as  they  are  discussed. 

Sequence  of  Events.  Air,  which  carries  minute  dust  particles,  is 
forced  into  the  precipitator  by  some  device  such  as  a  fan  or  by  the  forward 
motion  of  an  aircraft-mounted  version  of  EARC-I.  After  proceeding  some 
few  centimeters  into  the  entrance,  the  air  then  enters  the  charging 
section.  Here,  ionized  air  caused  by  a  corona  discharge,  imparts  a  negative 
charge  to  the  aerosols  present  in  the  air  stream.  The  particles,  during 
and  after  their  accumulation  of  charge,  are  drawn  in  accordance  with 
Coulomb's  law  to  grounded  collection  plates  which  act  as  anodes. 

Corona  Discharge 

It  will  be  shown  that  the  swiftness  with  which  a  particle  is  collected 
depends  initially  upon  the  mount  of  charge  it  acquires.  The  source  of  the 
charging  process  is  the  corona  discharge.  The  corona  mqy  be  defined  as 
that  electrical  state  which  exists  between  two  bodies,  the  potential 
between  which  is  sufficiently  below  the  potential  required  for  a  generalized 
breakdown  (Sparkover)  of  the  dielectric  between  them.  For  example,  when  a 
potential  difference  is  applied  between  two  parallel  plates  in  air,  there 
is  a  uniform  field  present.  Increasing  the  field  potential  is  possible 
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until  it  reaches  approximately  30  KV/cm  which  is  the  electrical  break¬ 
down  point  for  atmospheric  air. 

A  non-linear  electric  a  eld  is  necessary  for  a  corona  discharge  to 
occur.  Consider  two  electrodes:  a  point  for  the  cathode  and  a  plane  for 
the  anode.  Application  of  a  potential  across  the  electrodes  causes  a  field 
which  has  a  high  value  at  the  point  and  a  lower  value  at  the  plane.  By 
applying  the  proper  potential  difference,  electrical  breakdown  will  occur 
in  the  gas  close  to  the  point.  This  localized  breakdown  constitutes  the 
corona  discharge  which  is  necessary  to  the  charging  process.  A  further 
increase  in  potential  would  result  in  generalized  breakdown,  effectively 
short-circuiting  the  electrodes. 

Considering  the  electrodes  found  in  EARC-1,  it  can  be  seen  that  a 
potential  across  the  wires  and  plates  will  result  in  a  localized  electrical 
breakdown  in  the  air  surrounding  the  wire.  Since  the  wires  are  the  cathodes, 
it  can  be  visualized  that  there  are  two  distinct  electrical  zones  in  the 
space  betweei  the  electrodes.  Inmediately  surrounding  the  negative  wire, 
positive  ions  tend  to  congregate}  while  the  remaining  space,  and  by  far 
the  major  portion,  is  occupied  with  negative  ions  in  transit  towards  the 
grounded  collection  plate.  Corona  discharge  in  EAHC-I  occurs  for  total 
ionization  current  as  high  as  20  ma  at  a  charging  potential  of  23  KV 
(Ref  1:39) •  Exceeding  these  limits  usually  results  in  complete  breakdown 
(sparkcrver). 

Particle  Charging  Mechanisms 

Bombardme nfc  Charging.  In  bombardment  charging  the  ions  farmed  by 
the  corona  move  along  the  lines  of  fbrce  of  the  electrical  field  between 
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the  cathode  and  anode.  Upon  collision  with  an  aerosol  traveling  with  the 
air  stream,  charge  i s  transferred  to  the  aerosol.  As  the  particle  charge 
increases,  a  sufficiently  high  Coulombic  repulsion  to  other  ions  is  reached 
which  prevents  the  farther  transfer  of  charge  to  the  particle.  At  this 
point,  the  parti  cl p  is  said  to  hold  a  saturation  charge.  The  following 
assumptions  were  made  to  arrive  at  a  mathematical  expression  to  solve  for 
the  saturation  charge  acquired  by  a  particle  (Ref  16:1187): 

1.  The  aerosol  particles  are  spherical, 

2.  Particle  diameter  is  much  less  than  the  distance  between 
particles  in  the  airstream. 

3.  Tiie  immediate  region  surrounding  a  particle  has  a  uniform 
ion  concentration  and  electric  field. 

The  assumption  made  here  as  well  as  any  other  assumptions  made  in 
the  theoretical  development  of  EARC-I  will  be  discussed  later  in  the 
light  of  how  they  effect  the  experimental  results  obtained  in  this  study. 

The  expression  for  saturation  charge  acquired  by  an  aerosol  by 
bombardment  charging  is  (iief  7:35)- 

T  i  in1  'll  5r‘ 

|J  +  *■  l  J  e  (3.1) 

where  qQ  *  saturation  charge  on  particle  ( el ectrons) 

K  r  dielectric  constant  of  the  aerosol 
E  =  electric  field  in  the  charging  section  ( statvolts/cm) 
r  ^-particle  radius  (cm) 
e  *  Electronic  charge  (esu) 

It  can  be  seen  that  for  a  given  irateiiel,  the  saturation  charge  is 
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directly  proportional  to  the  radius  squared  and  the  electric  field. 

Diffusion  Charging.  Charging  by  this  method  is  the  result  of  the 
thermal  motion  of  ions.  Charge  is  imparted  to  the  aerosol  during  random 
collisions  vdth  ions.  Diffusion  charging  is  dependent  upon  the  particle 
size,  air  temperature  and  the  time  the  particle  spends  in  the  charging 
section.  The  expression  fcr  diffusion  charging  is  (Ref  7s ip-)  • 


r 


rkT 


LU 


I  -b 


rrUe^rnt 


kT 


(3.2) 


where  q  »  particle  charge  (electrons) 

U  *  rms  thermal  velocity  (cm/sec) 
n  *  ion  density  (ions/cnP) 
r  ”  particle  radius  (cm) 
e  »  electron  charge  (esu) 
k  -  Boltzmann  constant 
T  ”  temperature  (°K) 
t  ■  time  (sec) 

The  assumptions  made  &r  the  derivation  of  equation  (3.2)  are  the  same 
as  for  the  derivation  of  the  bombardment  charging  equation.  In  addition 
it  was  assumed: 

iu  All  ions  which  reach  the  particle  are  attached  to  it. 

Combining  Diffusion  and  Bombardmm t  Charging.  The  total  saturation 
charge  from  both  charging  mechanisms  considered  is  not  simply  the  addition 
of  the  charges  acquired  by  each.  The  differential  expression  combining 
both  mechanisms;  follows  (Ref  7:1;3): 
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where  A  ■*  11.16  Y)Y"  £■ 

B  -  1.61  x  10"6  n 
C  -  5.79  x  10_lii  Yi/r%£ 

D  -  15.7  x  10^  nr2 
j  -  5.56  X  io6/r 
n  *•  ion  density 
e  =  base  of  natural  logarithms 

Equation  (3.3)  is  a  non-linear  differential  equation  which  was  solved 
using  the  IBM  7090  Digital  Computer  (Ref  13:88).  A  discussion  of  the 
computer  solution  is  presented  in  Chapter  V  which  covers  theoretical 
predictions  of  EARC-I  performance. 


Collection  of  Charged  Particle b 

Although  EARC-I  is  a  two-stage  precipitator,  collection  of  charged 
aerosols  is  not  restricted  exclusively  to  the  collection  section.  As 
soon  as  a  given  particle  enters  the  charging  section  it  immediately  begins 
accumulating  negative  charge.  Simultaneously  it  feels  an  attractive  fcrce 
from  the  grounded  plates.  A  charged  particle  then,  is  under  the  influence 
of  two  simultaneous  forces  at  all  times  during  its  presence  in  the 
precipitator.  There  is  the  force  of  the  airstream  in  which  it  is  traveling. 
There  is  also  the  Coulombic  attraction  from  the  oppositely  charged 
collection  plates.  An  expression  for  evaluating  the  particle's  eomponait 
of  velocity  normal  to  the  airstream  flow  can  be  derived  from  two  basic 
principles: 


1.  Stoke' s  law  describing  the  motion  of  a  small  sphere  in  a 
viscous  medium 


£7777 


(3-U) 
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where  F  -  drag  on  particle  (dynes) 

n  -  viscosity  of  air  in  poise 

v  *  terminal  velocity  of  particle  in  ci^/sec 
o 

r  ■  particle  radius  in  cm 

2,  The  electrostatic  force  on  a  charged  particle 

f=  (3.5) 

where  E  *  field  strength  ( statvolts/cm) 
e  *  charge  on  the  electron  (esu) 
q  *  charge  on  particle  (electrons) 

Assuming  that  the  particle  attains  terminal  velocity  instantaneously} 
equating  equations  (3.U)  and  (3.5);  and  solving  for  terminal  velocity,  one 
gets 

=  (cmstmt)  %  £/r  (3.6) 

which  is  in  the  direction  determined  by  the  field  aid  the  charge  on 
the  particle.  For  JSARC-I,  airflow  is  parallel  to  the  collection  plates 
which  are  mounted  horizontally.  Therefore,*/,  will  be  in  a  vertical 
direction.  For  particles  with  radii  less  than  1.5  u,  there  is  an 
increased  probability  that  thqr  will  pass  between  air  molecules.  This 
would  result  in  a  hi$ier  velocity  than  is  indicated  strictly  on  the 
basis  of  charge  acquired.  The  Cunningham  correction  factor,  CUN, 
compensates  for  this  phenomenon 
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COIU  — 


+ 


r 


+- 


-0. 934x 

0.4  l  £ 


(3.7) 


The  corrected  drift  velocity,  v^,  is  the  product  of  the  particle's 
terminal  velocity  and  the  Cunningham  correction  factor 

v .  *  v  x  CUN  (3.8) 

a  o 

Typical  values  of  CUN  and  v^  are  diown  in  Table  3-1.  Although  Table  3-1 
includes  particles  with  radii  smaller  than  those  considered  in  this 
study,  it  is  interesting  to  note  that  for  very  small  particles,  the 
drift  velocities  are  greater  than  particles  with  larger  radii.  This 
is  contrary  to  what  one  would  expect  from  the  discussion  thus  far. 
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Table  3-1 


Typical  Values  of  Cunning am  Correction  Factors, 
Drift  Velocities  and  Charge  as  a  Function  of  Radius 


'Radius 

CUN 

Vw 

f  cm/sec') 

Charge 

(electrons) 

.12? 

2.02 

3.21 

16 

.20 

1.61 

2.80 

28 

.30 

1.39 

2.65 

46 

.425 

1.26 

3.09 

84 

.6 

1.19 

4.18 

169 

.85 

1.14 

5.65 

340 

1.25 

1.09 

7.95 

734 

2.0 

1.06 

12.4 

1880 

3.0 

l.Oli 

18.3 

4240 

4.25 

1.03 

25.5 

8490 

(Ref  7s72A) 
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IV.  Experimental  Procedtr  e 

The  basic  test  far  judging  how  closely  EARC-I  operates  to  the  theoretical 
model  was  a  comparison  of  the  size-distribution  of  particles  collected  during 
actual  operation  versus  the  sizes  predicted.  It  remained  to  collect  samples 
of  particulate  matter  actually  precipitated  from  the  atmosphere $  then  to 
make  a  particle  size-distribution  determination  on  those  samples.  This 
chapter  discusses  the  procedures  used. 

Collection  of  Samples 

There  are  many  methods  for  determining  the  size-distributions  of 
particles  in  a  sample  of  natter.  Examples  are  sieving,  sedimentation, 
elutriation  and  centrifuging  methods.  When  applied  to  the  investigation 
of  EARC-I,  each  of  the  methods  listed  has  some  objectionable  disadvantage. 

Each  of  the  indirect  methods  involves  first  collecting  particulate  matter  on 
the  collection  plates  then  devising  some  method  of  removal  from  the  plates 
fcr  the  actual  sizing  procedure.  It  was  reported  by  Baker  (Ref  1j14i)  that 
removal  of  collected  particles  is  a  difficult  problem.  At  best,  the 
intermediate  step  of  removing  precipitate  from  the  plates  wculd  involve  loss 
and/or  distortion  of  the  original  particles.  The  direct  sieving  method 
would  offer  the  least  objectionable  method  of  taking  samples  of  collected 
aerosols.  However,  sieving  would  necessarily  introduce  obstructions  to  air 
flow  or  perturbations  of  the  electrical  fieLds# 

Microscopic  analysis  3eemed  to  offer  the  best  method  of  examining 
collected  particles  as  they  were  actually  precipitated.  The  use  of  a. 
microscope  was  also  indicated  by  the  heterogeneity  of  the  particles 
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contained  in  atmospheric  air.  For  example,  a  sedimentation  process  would 
doubtlessly  dissolve  some  particles,  while  the  light  extinction  method  of 
size  analysis  would  be  invalid  for  the  wide  range  of  densities  found  in 
natural  aerosols. 

The  first  consideration  then  was  to  collect  the  particle  samples  in 
such  a  way  as  to  minimize  physical  and  chemical  changes  from  their  deposited 
state.  Accordingly,  it  was  decided  to  use  household  aluminum  foil  as  the 
collection  surface.  The  foil  was  tightly  stretched  across  the  surface  of 
the  collection  plates  and  taped  along  the  edges  to  minimize  geometrical 
deviation  from  a  plane  surface.  There  was  no  visible  adverse  effects 
apparent  during  operation  of  EARC-I  as  a  result  of  the  foil. 

The  collection  laboratory  (Figure  l;.l)  was  on  the  ground  floor  of  a 
two-stoiy  building.  Air  was  drawn  from  a  20  ft.  x  hO  ft.  open  courtyard 
in  the  center  of  the  building  and  was  exited  through  an  external  wall 
fifty  feet  away.  The  entrance  to  EAItC-I  is  a  six- foot-long  rectangular 
duct  which  passed  from  the  courtyard  through  the  laboratory  wall.  A  90° 
duct  elbow  is  installed  on  the  outside  end  of  the  duct  and  is  fitted  with  a 
1/2  in.  wire  screen  to  keep  out  large  objects.  A  gelman  thermal  anemometer 
is  mounted  in  the  exhaust  duct  to  measure  the  speed  of  the  air stream. 

Foils  were  exposed  to  the  precipitation  process  for  different  time 
intervals  in  order  to  get  a  sample  with  an  optimum  densl-ty  of  particles 
deposited  for  microscopic  viewing.  The  foil  exposed  for  one  hour  appeared 
to  give  the  best  density  for  that  purpose. 

After  a  precipitation  run,  samples  were  taken  immediately  following 
shut-down.  At  measured  distances  from  the  entrance  end  of  the  collection 
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Figure  U.l 

Arrangement  of  Collection  Laboratory 

1.  Exit  Duct  from  EAjKC-I 

2.  EAKC-I 

3.  Inlet  Duct  to  EAltC-I 

U.  Power  Supply  for  MRC-I 
5.  Courtyard 
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plate  (takai  as  station  0) ,  0.1?  mm- thick  cover-glasses  were  taped  to  the 
aluminum  foil.  The  microscope  samples  were  removed  by  cutting  around  the 
cover  glass,  removing  the  foil  and  coverglass  together.  Samples  were 
arbitrarily  taken  at  nine  collection  points.  A  collection  point  is 
defined  as  the  width  of  the  coverglass  minus  the  distance  of  the  coverglass 
along  its  edges  covered  by  the  tape  used  to  secure  it  to  the  foil.  Table 
U-l  shows  the  stations  included  by  each  collection  point. 

Microscopic  Analysis 

Since  the  camples  to  be  analyzed  were  collected  on  aluminum  foil,  it 
was  necessaiy  to  use  incident  limiting  for  viewing.  The  best  available 
instrument  for  this  type  of  work  was  the  Vickers  Projection  Microscope 
described  earlier. 

Magnification.  For  the  size  range  of  interest  in  this  study  it  was 
necessaiy  to  use  the  hipest  magnification  available  with  an  optical  system. 
This  implied  the  use  of  oil  immersion  lens.  The  objective  lens  used  throughout 
the  analysis  was  the  2.2  mn  focal  length  lens.  The  projection  eyepiece  was 
a  1$X  compensated  lens  with  the  mechanical  tube  length  set  at  230  mm  length. 

The  bellows  length,  which  determines  the  magnification  of  the  image  projected 
on  the  ground  glass  screen,  was  set  on  the  90  cm  position.  (See  Figure  2,3). 

The  lens  system  (which  is  takai  to  mean  the  object  lens,  projection 
lens,  mirror,  and  projection  screen)  was  calibrated  against  a  Bausch  and 
Lomb  object  micrometer.  The  micrometer  scale  was  projected  through  the 
lens  system  on  to  the  screen.  With  the  microscope  arranged  as  described 
above,  the  smallest  division  on  the  micrometer  known  to  be  0.01  ram, 
measured  $0.26  mm  +  l.00$  ran.  The  magnification  then  was  approximately 
5000X. 
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Table  li— 1 


Designation  of  Collection  Points 


Collect! on 
Point 

1 

2 

3 

U 

5 

6 

7 

8 
9 


Inclusive  Stations 
(cm) 

0.1  -  2.1 
5.5  -  7.5 
8.0  -  10.0 
12.0  -  lli.O 
15.0  -  17.0 
18.0  -  20.0 
21.5  -  23.5 
25.0  -  27.0 
30.0  -  32.5 
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Microscopic  Measurement.  There  are  several  methods  of  taking  the 
actual  measurement  of  particles  in  a  field  of  view.  Among  these  are 
msasurement  by  direct  observation,  projection  on  a  screen  and  photography. 
As  noted  above,  the  dickers  has  facility  of  projecting  the  image  on  a 
ground  glass  screen  for  greater  magnification.  This  is  sometimes  referred 
to  as  "empty  magnification",  but  for  the  purposes  of  sizing  in  the  range 
of  interest  here,  the  increased  magnification  is  validly  used.  The 
actual  measurement  consisted  of  measuring  the  Martin  statistical  diameter 
(to  be  defined)  in  all  measurements.  (See  Figure  U.2).  A  pair  of  dividers 
was  used  to  measure  the  projected  particle  diameter  on  the  screen.  The 
spread  of  the  dividers  was  then  measured  against  a  plastic  centimeter  ruler 
whose  smallest  division  was  a  millimeter.  This  method  is  analogous  to  the 
use  of  the  camera  lucida  for  measuring  particle  sizes.  Basically  the 
camera  lucida  is  a  way  of  projecting  a  magnified  image  onto  a  piece  of 
paper  where  it  is  traced  for  subsequent  measurement.  (Ref  2:?0),  Within 
human  error  measurements  were  read  to  the  nearest  0.25  mm  on  the  ruler  when 
measuring  the  spread  on  the  dividers. 

Small  Particle  statistics.  Statistics  as  applied  to  small  particles 
is  a  specialized  topic  which  is  most  thoroughly  treated  by  Herdan  (Ref  3). 

At  the  outset  of  a  microscopic  analysis,  one  must  reconcile  himself  to  the 
fact  that  he  is  literally  looking  at  a  "drop  in  the  ocean".  Nonetheless 
the  procedures  of  microscopy  are  longstanding  and  have  evolved  as  valid 
methods  for  representing  the  macroscopic  world  by  investigating  a  relatively 
very  small  sample. 

Two  of  the  more  common  statistical  diameters  used  when  measuring 
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Particle  Showing  Statistical  Diameters 
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particle  sizes  are  Martin's  diameter  and  Feret's  diameter  (Figure  In 2). 
Martin's  diameter  is  defined  as  the  chord  which  bisects  the  area  of  the 
image  presented  by  a  particle  viiai  viewed  through  a  microscope.  The 
bisecting  line  is  always  taken  parallel  to  a  fixed  direction  to  avoid 
introducing  bias  in  the  selection  of  any  measured  diameter.  Feret's 
diameter  is  the  mean  chord  perpendicular  to  two  opposite  tangents  to 
the  particle  outline.  The  tangents  are  also  drawn  in  some  arbitrarily 
chosen  direction.  In  this  study,  the  Martin  statistical  diameter  was 
measured,  taken  in  the  horizontal  direction.  It  has  been  shown  that  of 
the  two,  Martin's  diameter  has  less  error  than  Feret's  diameter  when 
considering  particle  surfaces  and  when  compared  to  the  Stoke' s  diameter 
(Ref  3s  1*6). 

The  next  logical  aspect  of  arriving  at  a  valid  statistical  size 
distribution  is  the  size  of  the  sasple  analyzed.  There  are  varying  opinions 
on  this  score.  DallaValle  indicates  that  for  most  measurements,  a  sample 
of  200  particles  is  sufficient  to  get  a  valid  representation  (Ref  2s69). 
Herdan  recommends  that  the  sample  contain  between  300  and  500  particles 
(Ref  3*1*7).  Table  IV  in  Skinner,  et  al  (Ref  12;9)  gives  the  minimum 
number  of  particles  to  count  for  a  probable  error  of  2%  as  1*00.  It  was 
therefore  decided  that  a  valid  representation  of  the  size  distributions  of 
the  samples  taken  from  EARC-I  would  result  from  counting  at  least  1*00 
particles  in  each  sample.  To  avoid  bias,  all  particles  in  a  given  field 
of  view  were  counted. 

An  important  assumption  which  has  been  illicit  in  the  discussion  of 
the  analysis  of  EARC-I  samples  should  be  stated  here.  It  is  assumed  that 
what  appears  as  a  particle  under  the  microscope  i.e.  a  unit  of  matter  whose 
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cohesive  forces  have  maintained  its  viewed  size  and  shape  is  and  was  an 
aerosol  originally  carried  in  the  airstream  being  drawn  through  the 
precipitator.  There  is  the  possibility  that  an  original  aerosol  upon 
entering  the  precipitator  could  have  been  altered  in  size,  shape  or 
chemical  consistency.  Such  changes  could  have  occurred  in  the  high  electric 
field  in  the  charging  section  or  upon  impact  with  the  collection  plate  when 
deposited.  There  is  also  the  possibility  that  what  appears  as  a  particle  is 
really  an  aggregate  of  many  smaller  particles. 

Designation  of  Particle  Groups.  It  would  be  meaningless  to  treat  each 
particle  measured  in  a  class  by  itself.  It  therefore  becomes  necessary  to 
group  arbitrarily,  particles  within  certain  size  limits.  Table  U-2  lists 
the  particle  groups  for  the  purposes  of  this  study  and  also  the  average 
radius  of  each  group.  It  is  necessary  to  use  an  average  radius  in  the 
computer  program  used  to  make  theoretical  calculations. 

Use  of  a  Particle  Size  Analyzer 

An  attempt  was  made  to  extend  the  range  of  particles  analyzed  during 
this  study.  For  this  purpose,  use  was  made  of  a  Particle  Size  Analyzer 
(PSA)  (Figure  U.3)  made  by  the  Southern  Research  Institute  for  the  U.  S. 
Public  Health  Service.  The  PSA  was  designed  to  measure  particle  size- 
distributions  in  atmospheric  aerosols  and  in  other  suspensions  of  solid 
particles  and  liquid  droplets  in  gases  (Ref  18:1).  It  covers  the  range  of 
particle  sizes  from  0.15  to  2.0  microns  diameter. 

EARC-I  was  modified  (Figure  i4.l1)  to  permit  monitoring  "grab  samples" 
of  the  airstream  with  PSA.  To  accomplish  this,  1/2  inch  holes  were  drilled 
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Table  l*-2 


Group 

Designation  of  Particle  Groups  by  Radii 

Inclusive  Radii  (u)  Average  Radius 

1 

1<-R*2 

1.5 

2 

2*  R<  3 

2.5 

3 

3£IUU 

3.5 

1* 

UR*5 

U.5 

5 

5^R<  6 

5.5 

6 

6SR<? 

6.5 

7 

7*  R<8 

7.5 

8 

8&R<  9 

8.5 

9 

9*R<  10 

9.5 
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Figure  U,3 

Particle  Size  Analyzer  (PSA) 

51  -  Selector  Switch 

52  -  On-Qi'f-Calibration  Switch 
R2  -  Eescrirainator  Potentiometer 

(Ref  18) 
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Figure  IwU 

il/UiC-I  Modified  for  Use  with  FSA 

1.  JiAitO-I 

2.  Plexiglass  Container 
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in  the  first  collection  plate  (figure  iuf>)  at  the  sane  collection  points 
at  which  samples  were  collected  previously  for  the  microscopic  analysis. 

A  plexiglass  container  (Figure  U.6)  was  constructed  to  capture  the  grab 
samples  tbr  monitoring.  By  making  a  size-distribution  determination  at 
successive  collection  points,  it  was  hoped  to  arrive  at  data  similar  to 
that  vhich  resulted  from  the  microscopic  analysis  perf armed  for  the  range 
of  larger  particles. 

The  attempt  was  unsuccessful,  hew  ever.  It  was  found  that  the  PSA 
was  adversely  effected  by  the  electric  fields  in  EARC-I.  As  the  power- 
supply  voltage  was  increased  from  zero  volts,  the  PSA  meter  began  to  give 
erratic  readings.  For  a  total  current  of  0.2  ma  and  above,  the  meter 
read  off-scale.  One  further  attempt  was  made  to  make  size-distribution 
measurements  of  the  air  immediately  before  its  entrance  into  EARC-I  and 
immediately  after  its  leaving.  Measur audits  at  both  locations  resulted  in 
meaningless  data.  It  appeared  that  the  high  velocity  of  the  air stream 
caused  the  PSA  to  give  erratic  readings.  Additionally,  at  the  exit  of 
EARC-I,  the  same  effect  as  noted  above  was  observed  when  the  applied 
voltage  was  increased  from  zero  volts. 

These  results,  although  unsuccessful,  are  included  as  possible 
information  for  other  investigators  in  this  area. 
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f 


2 


Figure  k.5 

1,  Collection  Plate  Modified  for  Use  with 
Particle  Size  Analyzer  (PSA) 

2.  Grid  from  Charging  Section  of  EARC-I 
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Figure  U.6 

Plexiglass  Container  for  Use  With  PSA 
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Figure  k.7 

!•  Plexiglass  Container  Mounted  on  EARC-I 
2.  Collection  Plate 
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V.  Theoretical  Predictions 

General 

The  variety  of  substances  of  which  aerosols  are  conqoosed  and  the 
wide  range  of  particle  sizes  present  in  atomspheric  air  have  made  it 
necessary  to  use  several  methods  simultaneously  to  investigate  the  size 
range  of  aerosols,  intensive  data  has  been  obtained  with  single  instruments 
from  which  it  is  possible  to  show  a  coherent  size  distribution.  Based 
upon  such  data,  Junge  has  shown  that  particles  in  the  range  of  0.1  to  10  u 
radius  have  a  size  distribution  which  can  be  closely  approximated  ty 

cLN _  *»  constant 

d(log  r)  rB  (5.1) 

where  N  ■  the  total  concentration  of  aerosol  particles 

(per  cra^)  from  the  lowest  size  limit  to  size  r 
B  =  a  value  between  2.5  and  3.5  depending  upon  geographic  location 
and  meteorological  conditions.  (Ref  5*9).  Taken  as  3.0  for 
theoretical  calculations  herein. 

It  is  assumed  for  the  purposes  of  this  study  that  equation  (5.1) 
does  in  fact  represent  the  size  distribution  of  natural  aerosols.  This 
assumption  makes  it  possible  to  compare  experimental  data  from  EAliC-I 
with  similar  theoretical  calculations. 

Computer  Analysis 

Stuart's  computer  solution  of  the  charging  equation  (3.3)  calculates 
the  charge  attained  by  a  particle  of  given  radius  at  any  time,  based 
upon  the  operating  parameters  of  EARC-I.  (Ref  13:12).  With  this  base,  a 
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program  was  written  (Appendix  A)  and  used  to  predict  where  each  particle 
of  a  given  group  would  be  deposited  by  the  precipitation  pi'ooess.  Junge'  s 
distribution  function  provides  the  theoretical  size  distribution  of 
particles  which  enter  the  precipitator. 

Use  is  made  of  Fortran  symbols  in  continuance  of  the  style  of  Stuart, 
which  makes  following  the  computor  program  easier.  Equation  (5*1)  is 
therefore  re-written  as 

d(EN)  -  A  .  1  d  (ER) 

HR**  RR  (S.2) 

where  EN  -  the  density  of  aerosols  (per  cm3)  from  the  lowest  size 
limit  to  the  radius  RR. 

A  ■  constant 

B  -  same  as  in  equation  (£.1). 

Letting  a  represent  the  lowest  size  limit, 

ti 

**  f  -9~> 

<t(M)  "A  (5.3) 

a, 

and 

EN  (RR)  -  EN  (a)  -  -A/_l  -  _1  ]  (5.W 

(RR3  A BJ 

where 

RR  =  radius  of  particle  of  interest 

EN(RR)  a  the  total  density  per  cm3  of  particles  with  radius 
a  to  pa. r tides  with  radius  RR. 

Fraction  of  Particles  Collected  ( FRTO) .  Since  the  particle  groups 
in  this  report  include  radii  in  a  range  of  values,  it  is  necessary  to 
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derive  an  expression  for  the  "density"  of  a  given  group.  The  density  of 
particles  in  Group  J  is  designated  OQNC(J),  -where  J  ■  1,  2,  ...9  kid  is 
called  the  concentration  of  group  J.  The  concentration  of  particles  in 
any  group  is  found  by  subtracting  the  density  from  equation  (5.U)  using 
the  smallest  particle  radius  in  that  group,  from  the  density,  using  the 
largest  particle  radius  in  that  group,  i.e. 


I  -fix)  **  *  /  -for)  *x  -  j  Ax) 

s-/  ft. 


(5.5) 


As  an  example,  consider  particle  Group  1  which  includes  particles 
with  radii  ranging  from  l.Ou  to  2.0u.  The  concentration  of  particles 
of  Group  (1)  is 

GONC(l)  -  Jen  (2)  -  EN  (a]]-  [EN  (1)  -  EN  (a)]  (5.6) 

which  becomes  from  equation  (5.U) 

GONG(l)  -  -A  jl  -  1  ~~1-  (~  A  V1  ~ 

B  [iB  RR(2)^J  V  [il 

and 

C0NC(1)  -  A  f  1  -  1  1 

B  [rR(1)B  RR(2)^J 

The  general  Fortran  statement  for  any  group  J  can  be  w 
CONC(J)  -  A  f  1  -  1  ~| 

B  |rr(j)b  HR(J+1)B  J  ^.8) 

where  CONC(J)  «•  the  number  of  particles  in  group  J  (per  cnr*). 

To  get  the  total  number  of  group  (J)  particles  which  pass  through 
EARC-I  during  a  given  sampling  run,  simply  multiply  the  group  concentration 
by  the  total  volume  of  air  sampled 

ODNC(J)  x  V  (5.9) 


HR(1)b J 


(5.7) 


itten 


35 


GNE/Phys/63-5 

where  V  -  total  volume  of  air  sampled  (cm^), 

It  fbllows  that  the  total  number  of  Group  (J)  particles  collected 
at  any  point  is  equal  to  the  total  number  which  passed  through  the 
precipitator  times  the  fraction  of  particles  collected,  JRTO 

FRTO(j)  x  CONC(J)  X  V  (5.10) 


where  SRTO(J)  -  fraction  of  particles  of  group  J  collected.  It  is 
important  to  note  how  fRTO(j)  was  derived.  The  time  is  computed  based 
on  the  particle's  drift  velocity.  If  the  particle  enters  the  precipitator 
this  close  or  closer  to  the  collection  plate,  it  will  be  collected.  If  it 
is  not  within  this  distance,  it  will  not  be  collected  and  will  leave  the 
precipitator  with  the  airstream.  (Ref  7:76). 

Addition  of  the  terms  implied  in  equation  (5.10)  gives  one  the  total 
number  of  all  particles  collected  in  groups  1  through  9 

HtTO(j)  x  CONC(J)  x  V  (5.11) 

J-l 

The  fraction  (of  th  ■  total  nuitber  of  all  particles  collected  from 
all  nine  groups)  which  is  in  any  group  J,  is  calculated  by  dividing 
equation  (5.10)  by  equation  (5.11): 


jRTO(j)  x  CONGO.  /  x  V 

9 

JRTO(J)  x  CONC(J)  x  V 

J-l 


(5.12) 


Equation  (5.12)  is  the  theoretical  equivalent  of  the  experimental 
data  shown  in  Chapter  VI.  It  is  noted  that  the  volume  of  air  sampled, 
V,  as  well  as  the  constant  A/B  from  equation  (5.7)  cancel  when  calcula- 
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ting  the  final  values  to  which  the  exp  eriment  data  is  conpared. 

Designation  of  Collection  Points*  There  were  a  total  of  nine  samples 
microscopically  analyzed  to  gather  the  experimental  data.  To  form  a 
comparison  for  theoretical  calculations,  the  following  assumptions  were 
made.  Since  each  slide  was  traversed  across  the  smaller  dimension,  i.e. 

2.0  cm,  a  collection  point  was  defined  to  be  2.0  era  wide  as  shown  in 

Table  it. 2.  For  the  purpose  of  the  computer  program,  the  fraction  collected 

at  any  given  station,  K,  is  the  fraction  collected  at  the  upper  limit  of 

the  station  minus  the  fraction  collection  at  the  lower  limit  of  the  collection 

point. 

For  example,  consider  the  fraction  of  particles  in  group  1,  collected 
at  station  5.  From  Table  k.2  one  sees  the  collection  point  5  extends  from 
station  15.0  cm  to  17.0  cm.  Therefore  for  this  exanple 
FRT0(1,5)  -  FRTOU,  17  cm)  -  PHT0(1,  15  cm). 

The  computer  output  is  MtT0(J,  K) 

Where  J  =  particle  group  number  1,2,3.... ...9 

K  =  collection  point  number  1,2, 3,.... 9 
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VI.  Results 

The  results  of  the  microscopic  experimental  investigation  and  the 
theoretical  digital  co input er  computations  are  shown  in  the  graphs  and 
tables  which  follow.  A  brief  explanation  precedes  each  group  of  data 
presented. 

Experimental 

The  experimental  data  represent  what  was  actually  observed  in  the 
microscopic  viewing  of  the  samples  collected  at  the  nine  collection  points. 
Table  6-1  and  figures  6.1  through  6.9  show  the  maber  of  particles  in  each 
group  observed  at  the  nine  collection  points.  Table  6-2  is  an  expression 
of  the  same  information  in  terms  of  the  percent  of  the  total  sample  considered 
at  each  collection  point  i.e.  the  UOO  or  more  particles  which  were  sized. 

For  example,  of  the  ljli6  particles  sized  at  collection  point  1,  there  were 
317  which  fell  into  group  1  (i.e.  317  particles  had  radii  between  1.0  u  and 
2.0  u).  That  is  71.0(3$  of  the  sample  considered  at  collection  point  1 
were  in  group  1.  Table  6-2  shows  relative  amounts  for  comparison  with 
theory. 

Qualitative  Observations.  During  an  experimental  investigation  there 
are  certain  observations  made  by  the  experimenter  which  are  not  amenable 
to  a  numerical  classification.  This  situation  results  from  such  con¬ 
siderations  as  equipment  limitations.  As  was  mentioned  previously,  the 
theoretical  limit  of  resolution  for  the  optical  microscope  is  approximately 
0.2  u.  Therefore  attempts  to  measure  particles  in  that  size  vicinity 
would  result  in  highly  questionable  size  distributions.  Yet  it  is 
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Table  6-1 


Number  of  Parti  ales  Experimentally 
Observed  at  Each  Collection  Point 


Collection 

Point 

Particle  . 
Group  * 

K9 

1 

BtfSl 

| 

3 

8.0 

to 

10.0 

(cm) 

1* 

12. ( 
to 
11*.  < 
(cm! 

5 

15.0 

to 

17.0 

(cm) 

6 

18.C 
to 
20. C 
(cm) 

7 

21.5 

to 

23.5 
(cm) 

8 

25.0 

to 

27.0 

(cm) 

9 

30.5 

to 

32.5 

(cm) 

1 

216 

273 

293 

285 

321* 

2 

92 

83 

85 

75 

81 

100 

87 

90 

63 

3 

15 

61* 

1*9 

31 

15 

32 

35 

21* 

9 

h 

12 

29 

32 

17 

6 

17 

16 

■ 

Bi 

3 

■ 

17 

16 

■ 

7 

2 

2 

2 

6 

3 

9 

m 

3 

■ 

3 

2 

■ 

3 

7 

0 

6 

2 

0 

— 

0 

■ 

0 

0 

■ 

8 

3 

■ 

8 

0 

2 

■ 

■ 

■ 

E 

9 

0 

B 

■ 

0 

D 

0 

0 

0 

H 

SAMPLE 

TOTAL 

1*1*6 

1*09 

1*16 

1*06 

1*05 

1  . 

1*13 

1*25 

1*07 

1*06 
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Figure  6,1 

Number  of  Particles  Experimentally 
Observed  in  Each  Group 

ko 


Figure  6,7 

Number  of  Particles  Exp  erimentally 
Observed  in  Each  Group 
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Table  6-2 


Percent  of  Total  Sample  Expeiimentally  Observed 
from  Each  Group  at  Each  Collection  Point 


Collection 

Point 

Particle 

Group 

1 

2 

3 

■ 

■ 

6 

1 

8 

9 

1 

71.08 

1:7.92 

51.92 

67.2b 

72.3b 

62.23 

65.65 

70.02 

2 

20.63 

20.29 

20.b3 

18.1*7 

20.00 

2b.  21 

20.  b7 

22.11 

3 

3.36 

15.65 

11.78 

7.63 

3.70 

7.75 

8.23 

5.90 

1.72 

b 

2.6? 

7.09 

7.69 

1*.18 

l.b8 

b.12 

3.76 

0.98 

0.98 

5 

0.90 

b.l6 

3.85 

1.72 

1.72 

0.b8 

1.18 

0.b9 

0.b9 

6 

0.67 

1.U5 

1.68 

0.73 

0.25 

0.73 

0.b7 

0.2b 

0.7b 

7 

0.00 

0.98 

0.1*8 

0.00 

0.00 

0.2b 

0.00 

0.00 

0.25 

8 

0.67 

1.95 

1.92 

0.00 

0.b9 

0.2b 

0.23 

0.2b 

0.25 

9 

0.00 

0.2b 

0.21* 

0.00 

0.00 

0.00 

0.00 

0.00 

0.25 

APPROXIMATE 

TOTALS 

100$ 

10$ 

100$ 

100$ 

100$ 

100$ 

10$ 

100$ 

100$ 

Example:  At  collection  point  1,  71.08  percent  of  the 
particles  sized  were  in  particle  group  1. 
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important  to  note  that  in  general,  the  number  of  particles  which  this 
observer  saw  below  1.0  u  in  radius  (i.e.  the  lower  limit  of  the  range 
of  interest  herein)  for  exceeds  the  number  of  larger  particles  recorded 
in  Table  6-1.  A  conservative  estimate  vould  be  to  state  that  there  were 
at  least  ten  times  as  many  particles  with  radii  less  than  1.0  u  as  there 
were  particles  counted  in  each  group. 

During  the  operation  of  EARC-I  sparking  was  observed  to  occur  in  the 
collection  section  as  well  as  the  charging  section.  The  sparking  in  the 
chirging  section  was  more  pronounced  than  that  in  the  collection  section 
which  is  estimted  to  be  at  least  a  power  of  ten  less  frequent  than  sparking 
in  the  charging  section.  This  observation  has  importance  when  considering 
the  phenomenon  of  "back  corona"  vhich  will  be  defined  and  discussed  in 
Chapter  VII. 

Theoretical 

The  theoretical  "data"  pres  on  ted  here  was  arrived  at  by  digital 
computation  (see  Appendix  A  for  program  listing).  The  theoretical 
equations  developed  by  Lamberson  for  particle  charging  and  fractions 
collected  were  manipulated  to  give  an  output  which  is  the  theoretical 
equivalent  to  experimental  data  presented  in  Table  6.2.  This  theoretical 
data  is  presented  in  Table  6.3.  The  operating  parameters  of  EARC-I  which 
were  input  for  the  computer  program  are  listed  below  in  Fortran  symbols: 


VL  - 

516.12 

cit/sec 

FRAR 

«=  1023.75  cm2 

SEP  ■= 

2.00 

cm 

E 

*  9750  volts/cra 

DL  - 

58.50 

cm 

CUR 

«  0,0012£>  anp 

DCH  = 

12.7 

cm 

AREA 

-  331.1:7  cm2 

SEPC  = 

2.0 

cm 

ECOL 

=  7600.  volts/cm 
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Table  6-3 


Percent  of  Total  Sample  Theoretically  Predicted 
from  Each  Group  at  Each  Collection  Point 


Collection 

Point 

Particle 

Group 

1 

2 

3 

■ 

5 

6 

7 

8 

9 

1 

78.50 

79.79 

81. 81* 

89.19 

100. 

100. 

0. 

0. 

2 

12.73 

12.92 

13.27 

10.80 

0. 

0. 

0. 

0. 

0. 

3 

U- 27 

i*.3i* 

1*.1*5 

n 

B 

B 

0. 

0. 

_ 

0. 

k 

1.91* 

1.97 

0.1*2 

a 

B 

B 

0. 

B 

0. 

$ 

1.0!* 

0.83 

n 

o. 

B 

0. 

B 

0. 

0. 

6 

0.62 

0.19 

o. 

B 

E 

0. 

B 

0. 

0. 

7 

0.1*0 

0. 

0. 

0. 

0. 

0. 

B 

B 

0. 

8 

0.27 

0. 

0. 

0. 

0. 

0. 

0. 

B 

0, 

9 

0.19 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

TOTALS 

100| 

100$ 

100$ 

100$ 

100$ 

100$ 

0$ 

0$ 

0$ 

Example:  Of  the  total  particles  collected  at 
collection  point  1,  78.30  percent  of 
them  should  be  in  particle  group  1. 
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Finally,  for  comparison  purposes,  a  log-plot  of  the  theoretical  and 
experimental  fractions  of  particles  oollected  in  each  group  at  each 
collection  point  is  presented  in  Figures  6.1  through  6.18. 
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Figure  6,11 

Experimeital  and  Theoretical 
Percent  of  Sarple  vs  Average 
Particle  Radius 
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2  5 

PARTICLE  RADIUS  (u) 


Figure  6.12 

Experimental  and  Theoretical 
Percent  of  Sample  vs  Average 
Particle  Radius 


QNE/Phys/63-5 


Figure  6,13 

Experimental  and  Theoretical 
Percent  of  Sample  vs  Average 
Particle  Radius 
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PARTICLE  RADIUS  (u) 


Figure  6.14 

Experimental  and  Theoretical 
Percent  of  Sample  vs  Average 
Particle  Radius 
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Figure  6.16 

Sxperiruental  and  Theoretical 
Percent  of  Sample  vs  Average 
Particle  Radius 
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PARTICLE  RADIUS  (u) 
Figure  6.17 


Experimental  and  Theoretical 
Percent  of  Sanple  vs  Average 
Particle  Radius 
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Figure  6,18 

Experimental  and  Theoretical 
Percent  oi’  Sample  vs  Average 
Particle  Radius 
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VII.  Analysis  and  Conclusi  ons 

General 

From  an  infection  of  Tables  6-2,  6-3  and  Figures  6.10  through  6.18 
one  can  see  that  there  are  two  extremes  of  agreement  between  the  curves 
representing  theoretical  predictions  and  those  representing  experimental 
observations.  On  the  one  hand,  there  is  exceptionally  close  agreement  far 
collection  points  1  through  J>.  That  is,  referring  to  the  appropriate  figures, 
one  can  see  that  for  the  first  5  collection  points,  the  maximum  deviation  of 
experimental  data  from  theory  are  as  shown  belowt 


Collection 

Point 

%  Deviation 
from  Theory 

Average  Radius  for  Which 
Maximum  Deviation  Occurs 

1 

56.0 

2.0 

2 

ItO.O 

1.5 

3 

36,6 

1.5 

it 

24.7 

1.5 

5 

27.0 

2.5 

where  %  Deviation  from  theory  _  theoretical  $  experimental#/ 

(for  given  radius  "  of  sample _  of  sample  /  (7.1) 

theoretical  %  of  sample 

In  general  it  can  be  seen  that  there  is  a  larger  deviation  from  theory  for 
particles  with  smaller  radii,  for  particles  in  the  size  range  considered  in 
this  study. 

On  the  other  hand,  beginning  with  collection  point  6,  there  is  an 
obviously  increasing  gap  betweai  the  experimental  and  theoretical  curves. 

At  collection  points  7,  8  and  9  there  is  an  undefined  deviation.  Where 


theory  predicts  that  no  particles  should  be  collected,  microscopic  analysis 
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shows  that  particles  have  been  collected  from  all  groups  of  interest. 

It  must  be  noted,  also,  that  the  curves  shown  in  Figure  6., 10  through 
6.18  are  relative.  They  do  not  show  a  count  of  particles  at  the  last 
collection  point  with  reference  to  the  first.  That  is,  one  cannot  say 
from  the  experimental  data  collected,  what  percent  of  total  number  of 
particles  which  alter  the  precipitator,  remain  in  the  airstream  by  the 
time  collection  point  9  is  reached.  As  an  illustration,  consider  a  test 
volume  of  air,  containing  N  particles  of  all  sizes,  entering  EAHG-I.  As 
that  volume  progresses  downstream,  particles  are  continually  being  removed 
by  electrostatic  precipitation.  As  the  volume  of  air  arrives  at  collection 
point  9,  for  instance,  all  that  can  be  said  about  the  test  volume  is  that 
it  contains  less  then  the  original  N  particles.  It  is  quite  possible  then  that 
at  the  later  collection  points  the  precipitation  is  working  on  a  veiy  small 
percent  of  the  original  number  of  aerosols  which  entered  with  the  airstream. 

Analysis  of  the  results  can  be  conveniently  discussed  in  two  parts. 

First  there  are  the  influences  which  occur  from  the  time  a  given  particle 
enters  the  precipitator  until  it  is  initially  deposited.  That  is,  until 
the  time  the  particle  first  contacts  the  co  llection  anode.  Secondly  th  ere 
are  phenomena  vhich  occur  after  a  particle  is  deposited  and  which  must  be 
considered  in  relation  to  the  experimsital  results  observed. 

Influences  Prior  to  Initi  al  Deposition 

Dielectric  Constant  (K).  In  the  equation  (3.1)  far  bombardment  charging, 
the  dielectric  constant,  K,  enters  the  expression  in  the  factor 


G  =  1  +  2(K-1) 

TTTT) 
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(Ref  7:35)«  Lamber son's  calculated  mean  value  of  K  ■  1*  by  averaging 
materials  usually  found  in  urban  areas  is  reasonable  for  arriving  at  a 
working  value  (Ref  7s 21).  However,  for  the  specific  case  under  consideration 
in  this  study,  one  must  allow  for  K  to  vary  between  a  value  close  to  zero 
and  very  large  values.  In  which  case,  G  can  be  seen  to  vary  between  0  and 
3. 


It  can  then  be  seen  that  the  drift  velocity  of  the  particle  is  directly 
proportional  to  G,  because  from  equation  (3.1) ,  saturation  charge 

q0  -  G  E,  r2  (7.3) 

The  drag  on  the  particle  from  Stoke’ s  Law  is 

F«6lTrnV0  (7.W 

and  the  electrostatic  force  between  the  particle  and  the  collection 
anode  is 


F-%eE?  (7.50 

Equating  the  expressions  for  drag  and  electrostatic  force,  and  substituting 
for  qQ,  the  drift  velocity  (without  Cunningham  correction)  is 

’  ‘  0  *1  V  (7.6) 

- 5 


for  a  given  set  of  precipitator  parameters.  It  can  be  seen  that  two 
particles  with  the  same  radius,  but  with  dielectric  constants  at  opposite 
ends  of  their  extreme  values,  could  be  deposited  over  a  relatively  wide 
range  of  collection  points  downstream  of  •the  entrance.  At  the  extreme 
values  of  X,  the  relative  collection  point  of  the  particle  with  the 
small  X  would  be  three  times  the  distance  to  the  collection  point  of  the 
particle  with  the  large  X. 
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It  is  interesting  to  note  further  that  as  the  particles  with  smaller 
radii  are  considered,  the  separation  of  collection  points  is  more  pro¬ 
nounced.  This  is  evident  by  looking  again  to  equation  (7.6)  and  noting 
that  the  drift  velocity  is  also  directly  proportional  to  the  radius. 

A  1.0  u  particle  deposits  farther  downstream  than  a  10.  u  particle,  for 
instance.  This  may  explain  in  part  the  larger  deviation  between  theory 
and  experimental  at  the  smaller  radii. 

Effects  of  the  field  Distribution  in  the  Charging  Section.  In  the 
derivation  of  the  equation  for  charging  an  aerosol  by  corona  discharge, 
it  was  assumed  that  the  electric  field  was  a  constant.  The  charging 
field  is  taken  to  be  the  value  of  the  voltage  potential  between  grid  wire 
and  ancde,  divided  by  their  separation.  This  is  really  the  expression 
for  a  linear  field.  A  corona  discharge  is  dependent,  however,  upon  the 
field  being  non-linear.  To  derive  an  analytical  expression  for  the  field 
potential  distribution  in  a  corona  Held  would  be  extremely  difficult  and 
has  not  been  reported  to  date.  Penney  and  Matick  used  a  probe  to  measure 
the  potentials  in  d.c,  corona  field  far  a  plate  to  wire  geometry.  The 
experimental  data  thus  obtained  substantiates  what  one  would  intuitively 
expect  the  field  to  be.  That  is,  the  field  would  be  greater  at  the  wire 
and  would  decrease  towards  the  plate.  A  characteristic  plot  of  equipo- 
tential  lines  are  shown  in  Figure  7,1. 

•It  can  be  seen  by  inspection  then  that  particles  which  enter  the 
precipitator  in  ary  given  plane  parallel  to  the  plate  electrode  will  be 
under  the  influence  of  different  values  of  electric  field.  In  the  extreme 
case,  a  particle  entering  the  charging  section  along  a  wire  electrode 
would  see  a  field  at  least  1.5  times  that  of  a  field  a  particle  which 
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entered  in  the  ssrie  plane,  but  midway  between  two  adjacent  wire  electrodes. 

The  drift  velocity  would  be  greater  than  for  particles  closer  to  the  corona 
wire.  This  effect  was  evident  on  the  collection  plates  of  EARC-I.  On 
the  charging  section  end  of  the  plate,  there  was  visible  to  the  naked  eye 
a  layer  of  deposited  dust  forming  an  outline  of  the  wire  electrodes. 

Field  Distortion  Between  Charging  and  Collection  Sections.  Figure  7.2 
shows  a  scale  drawing  of  EARC-I  with  an  exploded  view  of  the  physical  point 
of  separation  between  the  charging  and  collection  sections.  Note  that 
between  the  grid  in  the  charging  section  and  the  negative  plate  in  the 
collection  section,  there  is  a  one-inoh  plexiglass  insulator  resulting  in 
a  one-inch  gap  between  the  grid  frame  and  the  negatively  charged  plate.  The 
field  distribution  is  conceptually  indicated  in  the  blowup  of  Figure  7.1. 

One  can  visualize  the  perturbation  which  a  given  particle  "sees"  compared 
to  the  orderly  field  present  elsewhere.  This  perturbation  most  likely  results 
in  a  contribution  to  the  "spread"  of  points  at  station  and  beyond  (6,7,...). 

Effects  After  Initial  Deposition  (Re-Entrainment) 

The  precipitator  theory  which  has  been  developed  here  describes  particle 
motion  until  it  is  initially  deposited  on  the  collection  plate.  However, 
one  has  no  guarantee  that  a  particle  will  not  lose  its  charge  on  the  collection 
plate  and  then  be  swept  back  into  the  air  stream.  Such  particles  mi^it  be 
deposited  farther  downstream  or  carried  out  of  the  precipitator  with  the 
exiting  air.  However,  caused,  the  re-introduction  of  particles  into  the 
airstream  is  temed  re-entrainment.  It  is  important  to  note  that  a 
re-entrained  particle  is  not  necessarily  lost,  since  it  can  be  redeposited 
downstream  of  its  initial  point  of  deposit.  Baker  (Ref  1:U0)  used  a  Gelman 
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Figure  7.1 

Typical  Equipotential  Lines  in 
Plane-to-Wire  Geometiy  (Ref  10:95) 


Utt/ee  QATH0i>£~ 
(fS.SK/) 
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Figure  7.2 


Side  View  of  EARC-I  with  Full-Scale  Exploded  View  of  Air  Gap 
between  Grid  aid  Negatively  Charged  Plate  and  True  Locations 
of  Collection  Points  2-6. 
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Air  Sampler  to  examine  air  immediately  as  it  left  EARC-I,  He  found 
negligible  dust  traces  on  a  2.0  micron  Gelman  filter  paper,  following 
U-hour  sampling  runs.  This  would  indicate  that  re-entrainment  for  2- 
micron  and  larger  is  negligible  in  the  sense  that  such  particles  escape 
from  EARC-I.  However,  the  possibility  exists  that  ''internal"  re- entrainment 
followed  by  re-deposition  could  occur  in  EARC-I.  Such  a  phenomenon  could 
also  cause  particles  to  be  observed  where  theory  would  preclude  their 
collection.  This  is  especially  pertinent  to  larger  particles  since  they 
should  plate  out  before  smaller  particles  in  the  range  considered  in  this 
study.  Re-entrainment  may  explain  the  larger  particles  experimentally 
observed  at  collection  points  7,  8  and  9. 

Formation  of  a  Layer  of  Collected  Particles.  As  a  means  of  determining 
the  efficiency  of  a  precipitator,  Lamberson  derived  an  expression  for  the 
fraction  of  particles  which  would  be  collected  for  any  given  group,  ERTO 
(Ref  7:73).  This  fraction  is  based  on  the  simple  geometric  relation  between 
the  distance  t'e  particle  travels  towards  the  collection  plate  while  it  is 
in  th e  precipitator,  and  the  total  interelectrode  distance: 

distance  particle  travels  towards  collector 
5HT0  -  electrode  during  the  time  it  spends  in  precipitator  (7.7) 

distance  between  electrodes 

Considering  equation  (7.7)  in  relation  to  identical  particles  which  enter 

the  precipitator  at  the  same  distance  from  the  collection  plate,  it  must 

be  concluded  that  theoretically  each  such  particle  will  be  deposited  at 

the  same  point  downstream  of  the  entrance.  The  result  would  be  an  eventual 

"pile-up"  of  the  identical  particles  hypothetically  considered  here.  Such 

a  situation  in  EARC-I,  where  the  plates  are  used  dry  (not  treated  with  an 
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adhesive  to  retain  deposited  particles)  could  result  in  erosion,  by  the 
airstream,  of  the  particles  at  the  surface  of  the  collected  layer.  For 
particles  re-entrained  in  the  charging  section,  there  is  the  likelihood 
that  they  would  be  immediately  recharged  and  redeposited.  For  particles 
re-entrained  in  the  collection  section,  where  there  is  no  corona,  recharging 
is  not  possible.  The  particle  may  thai  leave  the  precipitator  with  the 
effluent  airstream  or  be  redeposited  if  it  has  retained  a  residual  charge. 

The  above  possibility  is  in  essence  a  case  of  re-entrainment  caused  by- 
physical  limitations  of  the  equipment.  That  is,  the  collection  and  retention 
of  an  unlimited  number  of  aerosols  at  a  given  point  is  not  a  logically 
acceptable  tenet. 

Back  Corona.  Another  cause  of  re-entrainment  is  the  phenomenon  known 
as  back  corona  or  reverse-ionization.  A  particle  which  has  been  deposited 
on  the  clean  collector  anode  will  give  up  its  negative  charge  and  assume 
the  polarity  of  the  collector.  The  result  is  a  repulsive  force.  The 
particle  will  be  repulsed  if  the  Coulombic  repulsion  is  greater  than 
the  molecular  forces  of  attraction  which  exist  between  the  particle  and 
collector.  For  very  small  particles,  molecular  forces  predominate  and 
retain  the  collected  particles  (Ref  8:540).  This  effect  was  noted  ty 
Baker  (Ref  1:1*4)  who  struggled  with  the  problem  of  removing  deposited 
particles  from  the  collection  plate  of  EARC-I,  in  order  to  concentrate  them 
into  a  sample  which  could  be  .used  for  radioactivity  analysis.  He  was 
successful  only  after  resorting  to  ultrasonic  cleaning.  Order  of  magnitude 
calculations  by  Lowe  (Ref  8:547)  indicate  that  larger  particles  with  high 
conductivity  will  be  repulsed  from  the  collector.  If  a  previously  deposited 

70 


GNE/Phys/63-S 

layer  of  particles  covers  the  anode,  a  voltage  potential  will  develop 
across  the  layer  far  resistivities  greater  than  lcA®  ohm-centimeters. 

Local  breakdown  occurs  (back  corona)  creating  ions  of  polarity  opposite 
to  the  cathode,  thus  neutralizing  the  original  action  of  the  precipitator. 
Since  the  entrance  to  EARC-I  is  fitted  with  a  wire  screen  to  keep  out 
extremely  large  pieces  of  airborne  debris,  the  back  corona  is  a  logical 
explanation  for  the  sparking  which  occurs  during  the  normal  operation 
of  EARC-I. 

Summary  of  Analysis 

It  is  evident  from  the  topics  discussed  above  that  there  are  many 
phenomena  which  can  enter  into  the  total  process  of  electrostatic  precip¬ 
itation.  The  factors  listed  above  are  listed  as  possible  explanations 
for  the  deviations  between  theoretical  and  experimental  data.  They  are 
offered  as  possiblr  explanations  only  because  an  investigation  of  the 
phenomena  themselves  was  not  an  object  of  this  study.  It  is  an  indication 
of  the  imperfection  of  the  model  introduced  by  Lamberson  (Ref  7)  that 
such  phenomena  cannot  be  treated  in  an  analytically  rigorous  manner. 
However,  in  spite  of  these  phenomena,  Lamberson1 s  simple  model  does 
give  some  correlation  as  shown  in  Figure  6,10  through  6.18. 

Conclusions 

The  experimentally  determined  size  distributions  of  precipitated 
aerosols  at  the  various  collection  points  when  compared  to  the  expected 
size  distribution  at  the  same  collection  points  lead  one  to  conclude 
that  EARC-I  does  collect  particles  roughly  in  accordance  with  Lamberson' s 
itcdel.  The  assumptions  made  in  the  original  derivation  of  the  theoretical 
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model  plus  the  variable  quantities  inherent  in  atmospheric  aerosols  not 
only  explain,  but  indeed  lead  one  to  expect  the  deviations  observed  in 
the  two  sets  of. data. 

The  two  basic  assumptions  are  here  restated  in  view  of  the  possible 
conclusions  to  be  drawn.  That  is,  the  theoretical  and  experimental  data 
should  agree  reasonably  wellj  if  the  size  distribution  of  naturally  occuring 
aerosols  is  that  described  by  Junge  and  if  EARC-I  does  operate  as  described 
by  Lamberson’s  model.  These  assumptions  are  reiterated  to  mention  that 
the  remote  possibility  exists  that  both  assumptions  may  be  wrong  and  that 
the  errors  in  one  cancel  the  errors  in  the  other.  This  possibility  is 
considered  to  be  unlikely. 
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■  VIII.  Recommendations 

In  spite  of  the  lack  of  a  complete  understanding  of  electrostatic 
precipitators,  their  use  has  increased  since  the  early  1900's*  Industry 
has  used  them  effectively  for  the  use  of  relatively  large  particulate 
matter  such  as  fly  ash  from  power  stations.  Baker  has  demonstrated 
conclusively  that  EARC-I  in  its  present  form  does  collect  more  than  adequate 
amounts  of  particulate  matter  for  analyzing  the  atmosphere  for  its  radio¬ 
activity  content.  This  is  true  although  the  method  he  used  for  secondary 
concentration  of  collected  samples  of  airborne  particulates  most  probably 
involved  much  loss  of  the  originally  precipitated  matter. 

The  literature  reviewed  by  this  author  during  the  course  of  this 
investigation  leads  him  to  conclude  that  the  theory  has  not  changed 
substantially  since  the  first  electrostatic  precipitators.  On  the  other 
hand  ingenious  engineering  has  rendered  them  useful  in  a  wide  variety 
of  situations,  in  fact,  even  where  theory  would  preclude  their  use 
(Ref  11:97).  It  is  therefore  recommended  that  the  development  of  EARC-I 
be  continued  from  an  engineering-development  point  of  view  rather  than 
to  pursue  further  theoretical  analyses. 

Observations  made  during  the  microscopic  phase  of  this  study  confirm 
the  thought  that  electrostatic  precipitators  are  especially  efficient  in 
the  precipitation  of  small  particles.  However,  for  the  purpose  of  radio¬ 
activity  detection,  one  cannot  consider  a  particle  properly  precipitated 
until  it  has  been  removed  from  the  atmosphere  and  conveniently  placed 
under  a  detector  for  analysis.  Therefore,  for  future  Air  Force  and/or 
Civil  Defense  use,  engineering  sophistication  of  EARC-I  should  be  effected 
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in  the  following  specific  areas,  listed  in  the  order  of  their  importance: 


1.  Find  or  devise  an  easily  removable  and  compressible  coating 
which  can  be  applied  to  the  collection  plates  of  EARC-I.  The  require¬ 
ment  for  compressibility  is  intended  to  include  any  means  of  making  a 
collected  sample  small  enough  to  be  used  with  existing  radiation  detectors. 

2.  Find  or  devise  a  method  of  making  the  coating  highly  conductive 
and  adhesive  in  order  to  render  the  coating  electrically  identical  to  the 
collection  plate.  This  would  lead  towards  100$  retention  of  deposited 
particles  because  of  the  combined  effects  of  the  Coulombic  attraction  and 
adhesive  quality. 

3.  Investigate  the  effects  of  reducing  the  collection  plate  area. 
The  ultimate  aim  is  to  attain  a  small  one-stage  precipitator. 

lu  Redesign  EARC-I  to  facilitate  the  removal  and  installation  of 
the  plates  and  grids. 

If  it  is  desired  to  continue  investigating  the  manner  in  which  EARC-1 
collects  microscopic  particles,  one  could  use  a  radioactivity  tracer  method 
in  the  following  general  manner: 

1.  Introduce  particles  of  known  size  and  material  (hence  constant 
K)  impregnated  with  a  known  radioactivity  into  the  inlet  air  stream. 

2.  Measure  the  level  of  radioactivity  at  various  points  along 
the  collection  plate. 
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*  Physical  Research  Laboratory,  28Q  Building,  Dow  Chemical  Company,  Midland, 
Michigan  has  available  Polystyrene  spheres  which  come  in  the  following 
sizes:  0.09  u,  0.19  u,  0.26  u,  0.37  u,  0.56  u,  0.8l  u,  1.1  u  and  3.0  u. 

It  also  has  the  following  size  mixtures  available:  6.0-lU.Ou,  16.0-28.0  u, 
29.0-56.0  u,  and  1*9.0-103.0  u.  Such  a  spectrum  of  sizes  would  allow  for  a 
wide  range  of  investigation. 
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3.  Calculate  the  number  of  particlee  required  to  give  the  level 
of  activity  detected  at  each  test  point. 

It.  Repeat  the  procedure  on  clean  plates  using  particles  of  the 
same  material  as  in  1  but  with  a  different  known  radius. 

The  procedure  suggested  above  would  correct  the  limitations  inherent 
in  sampling  naturally  occurlng  aerosols  as  was  done  in  this  investigation. 
Specif ically,  one  could: 

1.  Know  the  absolute  number  of  particles  introduced  into  the 
precipitator. 

2.  Know  the  dielectric  constant,  K 

3.  More  closely  approach  the  assumption  that  all  particles  are 

spherical. 

A  tracer  technique  would  also  lend  itself  to  investigating  re-entrainment. 
Re-entrainment  could  be  studied,  for  example,  by  pulsing  the  radioactive 
particles  into  the  entrance  and  observing  the  radiation  response  as  a  function 
of  both  distance  and  time  after  pulse. 


GNE/Fhy s/63-5 


Bibliography 


1.  Baker,  J.  W.  Sampling  of  Airbo rne  Radioactive  Particles  by 
Electrostatic  Precipitation.  M.S.  Thesis,  Air  Force  Institute 
of  Technology,  19 62. 

2.  DallaValle,  J.M.  Micromeritics  (Second  Edition).  New  York:  Pitman 
Publishing  Corporation,  1948. 

3.  Herdan,  G.  Small  Particle  Stati sties  (Second  Edition).  New  York: 
Academic  Press  Inc,,  1966. 

1*.  Hewitt,  G,  W.  llThe  Charging  of  Small  Particles  for  Electrostatic 
Precipitation. 11  Transaction:  American  Institute  of  Electrical 
Engines,  Communications  and  Electronics  -  Part  I,  75:  300-306 

5.  Junge,  C,  E.  "Atmospheric  Chemistry."  Advances  In  Geophysics,  4 
New  York:  Academic  Press  Inc.,  1955. 

6.  Junge,  C.  E.  "The  Chemical  Composition  of  Atmospheric  Aerosols  I: 
Measurements  at  Round  Hill  Field  Station,  June-July  1953." 

J ournal  of  Meteorology,  11:  323-333  (August  1954). 

7.  Lamberson,  D.  L,  A  Study  of  the  Electrostatic  Precipitation  of 
Radioactive  Aerosols.  M.~~S.  'Thesis,  Air  Force  Institute,  of 
Technology,  1961. 

8.  Lowe,  H.  J.  and  Lucas,  D,  H.  "The  Physics  of  Electrostatic 
Precipitation."  British  Journal  of  Applied  Physics,  Supplement  2: 
SUO  -  SU7  (1953). 

9.  Murphy,  H.  T.  jst  al.  "A  Theoretical  Analysis  of  the  Effects  of 
an  Electric  Field  on  the  Charging  of  Fine  Particles."  Trans¬ 
actions:  American  Institute  of  Electrical  Engineers ,  Communi¬ 
cations  and  Electronics  -  Part  I,  78:  318  -  325  (September  195 9) . 

10.  Penney,  G.W.,  and  Mat.ick,  R.  E.  "Potentials  in  D-C  Corona  Field; 
Transactions:  American  Institute  of  Electrical  Engineers,  Part 
79:  91~99TMay  i960). 

11.  Rose,  H.E.,  and  Wood,  A.J.  An  Introduction  to  Electrostatic 
Precipitation  In  Theory  and  Practice.  London:  Constable  and 
Company  Ltd. ,  1 9561 

12.  Skinner,  D.G.,  et  al.  Determination  of  Particle  Si?e  in  Sub- 
Sieve  Range.  London:  The  British  Colliery  Owners  Research 
Associ  at  ion,  1944 


76 


GW E/Phy s/63-5 

13»  Stuart,  R.S.  The  Analysis,  Optimization,  and  liesign  of  an  Airborne 
Radiation  Sampling  Electrostatic  Precipitator.  M.  S,  Thesis,  Air 
Force  Institute  of  Technology,  19t>2. 

lit.  White,  H.J.  "Modern  Electric  Precipitation,"  Indu strial  and 
Engineering  Chemistry  1;?:  932-939  (May  1555), 

15.  White,  H.J.  and  Penney,  G.W.  "Electrical  Precipitation  Fundamentals." 
Engineering  Proceedings  P-39  (July  1961) . 

16.  White,  H.J.  "Particle  Charging  in  Electrostatic  Precipitation." 
Transactions:  American  Institute  of  Electrical  Engineers  -  Part  II,  70: 
11U0-1191  U951TT 

17 .  No  Author.  1-hotomlcrot-raphy  with  the  g icker  Projection  Microscope. 

York,  England:  Cooke,  Trough ton~  &  Simms.  " 

18.  Ho  Author.  "Reference  Manual  for  the  Southern  Research  Institute 
Particle  Size  Analyzer"  Birmingham. 


77 


GNK/rhys/63-5 


Appendix  A 

Fortran  Program  Used  for  Theoretical  Computations 


General 

The  IBM  7090  Digital  Computer  was  used  to  solve  the  theoretical 
equation  describing  the  motion  of  particles  as  they  flow  through  EARC-I. 

A  description  and  listing  of  the  Fortran  program  used  follows.  The 
dictionary  of  variables  is  listed  elsewhere  because  the  text  of  this 
report  uses  Fortran  symbols  where  possible  to  facilitate  understanding 
the  program.  Since  this  is  a  modified  version  of  Stuart's  program  (Ref  13) 
which  is  explained  thoroughly  in  his  thesis,  only  a  brief  program  description 
is  offered  below. 

Program  Description 

The  purpose  of  this  program  is  to  compute  the  fraction  of  particles  in 
each  group  collected  at  each  collection  point  based  upon  the  theoretical 
equations  described  in  previous  chapters.  The  program  also  has  the 
convenience  of  computing  experimental  fractions  so  collected,  based  upon 
the  input  of  experimental  data. 

The  basis  of  the  program  is  found  in  statements  number  691+1  through 
statement  number  35.  This  section  of  the  program  computes  the  total 
charge  a  particle  has  attained  at  any  time  from  t  ■  0  seconds  to 
t  3  1001  seconds  at  0.0001  -  second  intervals.  From  the  charge  QQ  the 
drift  velocity  of  the  particles  of  interest  is  computed  for  ary  timej 
and  as  described  in  the  text,  the  point  at  which  the  particle  is  finally 
theoretically  deposited.  Statements  number  799  through  80U  accomplish 
the  latter  computations.  The  program  is  then  cycled  nine  times,  once 
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far  each  particle  group  then  eight  times  on  top  of  that.  The  latter 
ei $it  cycles  are  arbitrary.  Thqy  are  for  the  purpose  of  observing  the 
effects  by  varying  the  exponent  BETA  of  Junge's  differential  equation 
for  aerosol  distribution  in  the  atmosphere. 

The  output  of  the  program  is  a  listing  of  the  theoretical  prediction 
of  the  fraction  of  each  particle  group  collected  at  each  collection  point 
for  each  value  of  BETA  investigated.  The  output  is  in  the  form  of  an 
arrsy  the  rows  representing  particle  groups  1  through  9;  the  columns 
representing  collection  points  1  through  9.  All  output  is  labeled  with 
an  appropriate  heading. 

Fortran  Program  Listing 


Hff  AM/uble  copy 
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10  COMPUTAT I ON  OF  FRTO(J,K)  USING  AVER  AC-F  P.ADIU1.5,  2.5.  ...9,  5) 

DIMENSION  G  i  1001  )  ,  R  (  9  )  ,  C  U  N  (  9  )  f  V  (  1 0  0  L  )  »  C  CiN  C  (  9  ) 

DIMENSION  TSTNf  18)  ,STN(  IS)  ,  FRT0<9, 19)  ,C0LM(91  ,EN(  10)  ,KR(  10) 
DIMENSION  F  R  (  9  )  ,  F  R  C  (  9  »  9  )  , SEN  ID)  ,  E  F  F  (  9  ) 

DIMENSION  EXPERI9,9),FXC0L(9),XFRT0(9t9i 
100  FORMAT (Ell. A) 

115  FORMAT  (5X, 3 'INFRACTIONS  COLLECTED  -  THE  OP.E  T  I  C  A{_ ,  /  /  t  9E  1  4 . 6/ )  ) 

120  F0RMAT(5X,7HBETA  =  E14.3/! 

130  FORMAT ( 5X, 3oHFRACT IONS  COLLECTED  -  EXPER I  MENTALLY ,//( 9E L4. 5/ )  ) 

131  FORMAT (5X»6HCUN<J) » / / ( 9E l 4 . 6/  !  ) 

132  FORMAT ( 5  X  ,  4  8HNUMB E R  OF  PARTICLES  IN  EACH  GROUP  -  EXPERIMENTAL ,//( 9 
IE  14.6/  )  ) 

133  FORMAT I5X*  24HM.AS  S  FRACTIONS  PER  GROUP , / / ( 9E  l  4 . 6/ )  ) 

2  READ  INPUT  1APF  2 , 1 00 , ( R ( J ) , J  =  1 , 9 ) 

3  READ  INPUT  TAPE  2 , 1 00 , VL , SEP , DL . DCH , SE PC , FR AR , E , CUR , ARE  A , EC OL 
500  READ  INPUT  TAPE  2 , 1 00 , ( S TN ( K ) , K= l ,  l  8 ) , ( RR ( K  )  , K =  1 , 1 0 > 

READ  INPUT  TAPE  2 , l 00 , ( ( E X P ER ( J , K ) , J = U 9 ) , K = 1 f 9 ) 

5000  DO  5001  J  =  1  »  9 

5  00  1  CUN(J)=l.  +  ((9.42E-6)*((l.?3)  +  (.41*CXPF<(-0.93  4)*(10.*M5.»R(J))))) 
L)  )/(R< J)  )  | 

DO  5002  K-1,9 

EXCOL ( K ) =0.  j 

DO  5002  J  =  1 , 9 

500  2  EXCOL  (K)  =  EXCOL(K)t EX PCR! J,<) 

DO  5003  J  =  l,9 
DO  50C3  K  =  I  »  9 

5  00  3  XFRTOI J,K)=6XP5R( J , K ) / EXCGL  (  K ) 

DO  5005  J=  1,9  ! 

5005  FR  (  J  )  =LOGF  (  RR  (  J 1  )  /RR  (  J  )  ) 

WRITE  OUTPUT  TAPE  3 , l 3 1 , \ CUN ( J  )  , j =  l  , 9 ) 

WRITE  OUTPUT  TAPE  3 , 1 32 , ( ( FXPER ( J , K ) , X  = 1 , 9 ) , J= 1 , 9 ) 

WRITE  OUTPUT  TAPE  3 , 1 3 3 , ( FR ( J  )  , J = 1 , 9 ) 

BETA=1. 

DO  2  000  N=  I  ,8 

WRITE  OUTPUT  TAPE  3, 120, BETA 
DO  69  K= l , 10,  l  I 

6  9  E  N  (  K  )  =  1 .  /  R  R  (  K  )  *  #  B  F  T  A 

DO  691  J= 1 , 9  ! 

691  CONC ( J ) =FN( J ) -EN( J  +  L  ) 

J  =  1 

r>ELT  =  0.0001  | 

VLCL-=  VL 

G  =  (  .  3  5  2  E  1  9  *  C  U  R  )  /  (  E  *  A  R  E  A  ) 

4  H  =  G*R( J ) **2  .  *E 
A= 1 1 . 16*H 

B  =  2.*A*  7.203E-8/ ( E*R<  J )**2.  } 

C-B*7.203E-B/(2.*C*RIJ)**2.) 

0=15.7E4*G*R(J}**2. 

F=5. 56E-6/R ( J ) 

T=DCH/VL 

I T= (  ( T/DELT  H-1.0) 

I  =  l 

Q ( I ) =0. 

5  X YZ=F*G ( I ) 

IFIXY  7.  -  87.4)9,9,7 
7  QP= A-3*Q ( l ) +C*Q ( I )*Q( I ) 
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10  computation  of  f  r  r  o  (  j  ,  k  i  using  average  rapiic.5,  2.5 - 9.5) 

GO  TO  L  5 

)  I F { XYZ-. 00001 >11,11,13 
1  l  Q  f’  =  A  -  P.  »  0  (  I )+C*G( I )*Q([)+D 
GO  TO  15 

1  3  QP=A-3»Q(  I  1 +C  *Q (  1  )  *G(  I  )+D*EXPF(-XYZ  ) 

15  1  =  1  +  1 

IF( 1-1001)1 7,17, 20 
170(1 ) =Q ( I-l )+QP*DELT 
GO  TO  5 
20  T  -0  . 

2  5  T  =  [ )  C  H  /  V  L 

I =( T/UELT  )  +  l  . 

Y  =  G (  I  +  l  ) -0 (  I  ) 

U  =  1  -  l 

U  =  1J*DFLT 
X  =  T  -  U 

ni<  =  Y«  X/OELT 

35  00=0  ( I )  +  n i ’> 

709  DO  POO  K=  1  ,  1  H 
793  OAC-O. 

TSTN(K)=ST‘J(K)/VL 
TnCH  =  ('CU/VL 

I TSTN  = ( TSTNI K ) /DEL  T  )  +  l  . 

I T  DCH= ( TDCH/DELT )  +  1. 

IF(STN(K)-DCH )33,33,34 

3 3  DO  8099  I  =  1 , l T  5  TN , 2 

V(  I  )=1.602E-12*C(  I  )*E*CU N(J)/(6. 0*3. 1416*1. 3E-04*R(J)  1 
3000  DAC=DAC+ (VI  I  1 *  D  E  L  T  1 

3  3  V  rCR= l.60?E-l?*CQ*E*CUN( J 1/(6. 0*3. 1416*1. SE-04*R<J1) 

999  0  A  C  =  D  A  C  + ( V  T  E R  *  X  1 

FCOOAC/SLPC 
IFIFCC-1; 1300,3003,3003 

34  00  SC01  1  =  1,1  TDCH »  2 

V  (  1 1 =1 .602C-I ?*G(  l 1 *E»CUN( J 1 / ( 6.0»3. 14  16* 1. 3E-04*R( J } ) 

800  1  DAC=DAC+  (  V(  I  1 »Df L  V 1 

VTFR  =1.602F-12»QG»E*CUN(J)/(6.0*3.1416*1.3E-04*R(J1) 

DAODAO  (  V T C R * X  1 
F  C  C  =  D  A  C  /  S  t  P  C 
IFIFCC-1.  1  122,3003,3003 
122  Vll  =  4. 73F  -  10*G0*ECQL*CUN(  J  )/R(  J  1 
D  A  = ( V  D 1 » !  ( STOCK l-DCHl/VL) 

F  C -  r )  A / ( SEP-OAC ) 

F  C  C  =  F  C  C  + (  (  l .-FCC 1 *FC  ) 

IFIFCC-1. 1300,3003,8003 
900  3  F CC  = 1 . 

300  F  R  TO ( J , K 1 =  FCC 
J  =  J  +  1 

IF(J-9) 4,4,801 
801  DO  802  J=l,c> 

DO  302  K=l,9 

80  2  FR  TO ( J , K 1 =FR TO < J , 2*K 1 -FRTOI J , 2»K- L 1 
00  8  06  J  =  l  ,  ) 

DO  9  0  6  K  = l , 9 

89  6  F  R  T 0  (  J  ,  K  1  =  F  ",  T 0  (  J  ,  K  )  «CONC  (  -J  1 
DO  80  3  K=l,f> 

COLY(K) =0. 
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10  COMPUTATION  OF  FRTOIJ.K)  USING  AVERAGE  RADIK1.5,  2.5. 
1)0  803  J=1,<J 

803  C0LM(K]=C0LM(K)4-FRT0(  J,K) 

DO  904  J=  L  » *) 

DO  804  K=L,9 

804  FRTOI J,K)=FRTO{ J,K}/COLM«K) 

805  WRITE  OUTPUT  TAPE  3 , 1 15 , U FRTO ( J , K ) ,K= 1 , 9 ) , J= 1 , 9) 

2000  3ETA=8ETA+.5 

5004  WRITE  OUTPUT  TAPE  3 , 1 30 , I ( XFRTOI J , K ) , K= 1 , 9 ) , J= l t 9 ) 

CALL  EXIT 

ENDI1 ,0,0, 0,0 .0,1 i 0.0, l ,0.0. 0.0.0) 
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